Novel vortex structures in dipolar condensates 
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We investigate the properties of single vortices and of vortex lattice in a rotating dipolar conden- 
sate. We show that vortices in this system possess many novel features induced by the long-range 
anisotropic dipolar interaction between particles. For example, when the dipoles are polarized along 
the rotation axis, vortices may display a crater-like structure; when dipoles are polarized orthogonal 
to the rotation axis, vortex cores takes an elliptical shape and the vortex lattice no longer possesses 
hexagonal symmetry. 
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Quantized vortices represent one of the fundamental 
topological structures in nature. Studies of quantized 
vortices were launched out of the efforts to understand 
the flow properties of superfluid 4 He, but have now ex- 
tended into a variety of systems including superconduc- 
tors, superfluid 3 He, quantum magnets, liquid crystals, 
nuclear matter and even cosmic strings 0. The realiza- 
tion of weakly-interacting quantum gases of atomic va- 
pors provides a new platform for the study of quantum 
vortices. Vortices in atomic Bose-Einstein condensates 
have been produced, perturbed and engineered by pre- 
cise external fields, and can be studied theoretically in 
a quantitative manner using first-principles calculations 
@ . Very recently, vortices were also observed in a quan- 
tum degenerate Fermi gas Q. 

So far, vortices have only been produced in conden- 
sates featuring short-range s-wave collisions. The cre- 
ation of condensation in the transition metal 52 Cr of- 
fers a completely new system of quantum gases Q. The 
ground state of 52 Cr has a large magnetic dipole mo- 
ment of \x = 6/ib, with [ib being the Bohr magneton. 
Hence the magnetic dipolar interaction between atoms 
plays a crucial role in 52 Cr condensates. Previous theo- 
retical studies of dipolar quantum gases have focused on 
the ground state properties of non-rotating systems p| . 
The purpose of this paper is to investigate the properties 
of vortices in a rotating dipolar condensate. As we will 
show, due to the long-range and anisotropic nature of the 
dipolar interaction, single vortices and vortex lattices in 
such systems possess many unique features not present 
in condensates with only contact collisions. 

Model — We consider a trapped condensate of N 52 Cr 
atoms whose dipole moments are polarized along certain 
direction. The two-body dipolar interaction potential is 
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where fiQ is the vacuum magnetic permeability, and 9 
is the angle between the dipole moment and the vector 
r — r'. The trapping potential is assume to be harmonic 
with axial symmetry 



with u>± being the radial trap frequency and A = uj z /uj± 
the trap aspect ratio. For the sake of simplicity, we 
will focus on a pancake-shaped system with A ^> 1 
such that the condensate can be regarded as quasi two- 
dimensional (2D) whose motion along the z-axis is frozen 
to the ground state of the axial harmonic oscillator. 
Therefore the condensate wave function can be decom- 
posed as <t(r, t) — ip(p,t)(j)(z), where p = {x,y) and 
4>{z) = (A/V) 1 / 4 e~ Az / 2 . Throughout this paper, we 
adopt a unit system where the units for length, frequency 
and energy are given by a± = \J h/McJ±, and Huj±, 
respectively. The reduction of the effective spatial dimen- 
sionality considerably reduces the computational inten- 
sity. Although intrinsic three-dimensional (3D) phenom- 
ena such as the bending and Kelvin mode excitation of 
vortex lines cannot be studied, the 2D model does cap- 
ture many essential features of vortices at equilibrium, 
for example, the structure of the vortex core and that of 
vortex lattice. 

We assume that the condensate is rotating along the 
symmetry axis z with angular frequency f2. After in- 
tegrating out the z-variable, we obtain the effective 2D 
Gross-Pitaevskii equation in the rotating frame as 
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where fi is the chemical potential, L z the z-component 
of the orbital angular momentum operator, g^D — 
2(27rA) 1//2 7Va/aj_ the rescaled 2D collisional interaction 
strength with a being the s-wave scattering length which 
can be tuned via Feshbach resonance 0. The rescaled 
2D dipolar interaction potential takes the form 



VlF(P,P') = J dz\cf>(ztV dd (r,r') 



(4) 



U(r) = hlLul(x 2 +if + \ 2 z 2 ) 



(2) 



In the following, we will treat separately the cases where 
the dipoles are polarized along the axial and along the 
transverse direction, as they lead to very different phe- 
nomena. 
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FIG. 1: (Color online). Solid lines: function [/!(£>)]" 1/3 for 
A = 11, 100, and 400 in descending order. The dashed line 



represents the asymptote (2tt) 



-1/6, 



Axially polarized dipoles — For this case, we can find 
from Eq. J3J that 



where g = \p — p'\ and 

A 3/2 e A e 2 /4 



(5) 
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with i£n(z) being the nth-order modified Bessel func- 
tion of the second kind and U(a, 6, z) being the conflu- 
ent hypergeometric function. Fig. ^ shows the function 
[fx(g)]^ 1 ^ 3 for various values of the trap aspect ratio A. 
It can be seen that, as expected, the effective dipolar in- 
teraction is repulsive and isotropic. Independent of A, 
fx(o) asymptotically approaches ^/2~n q~ 3 at large g. At 
the other limit, as g — > 0, /f (f?) diverges much slower 
than g~ 3 , although the detailed behavior at small g de- 
pends on the value of A. Note that for a true 2D sys- 
tem with dipoles located in the xy plane and polarized 
along z axis, one expects a dipolar interaction potential 
Vj® cx g~ 3 as can be easily seen from Eq. (Q. For the 
quasi-2D case considered here, the effective dipolar inter- 
action potential deviates from Vj^ at small values of g 
in a fashion that makes the singularity of Vj® at g = 
integrable, a quite important property for numerical cal- 
culations. An alternative, and often more efficient, way 
to treat the dipolar terms in numerical calculations is to 
use the Fourier transform of V%P . 

The ground state wave function can be obtained by 
evolving Eq. in imaginary time. In the numerical 
results presented in the paper, we use N = 10 5 and u>± = 
2tt x 100Hz. First we want to study the structure of a 
single axial vortex at the center of the cloud. Due to 
the azimuthal symmetry, such a state can be written as 
ip(p) — x(p) ettp with ip being the azimuthal angle and x a 
real function satisfying the following ID radial equation: 
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FIG. 2: (Color online), (a) Radial density profile of the single 
vortex state for A = 100 and various scattering lengths (in 
units of Bohr radius, as), (b) The A dependence of critical 
scattering length (in units of as). Here □ and A denote, 
respectively, the non-vortex state and the axial vortex state. 



where 

/■OO />27T 

V(p)= p'dp' d^'V d 2 d D ( P 2 + P 12 - 2pp'cos V ')x 2 (p'). 
Jo Jo 

Figure Ufa) illustrates the radial density profile x 2 (p) 
of the single vortex state for different values of the scat- 
tering length. An interesting feature one can notice is 
that the density close to the vortex core oscillates as 
long as a is not too large. Such oscillatory behavior is 
induced by the dipolar interaction and gives the vortex 
a crater-like shape. Similar density oscillations are also 
observed in numerical studies with other forms of non- 
local interaction potentials which are originally employed 
to model the inter-particle interactions in superfluid 4 He 
0, . Note that no oscillations are found in the ground 
state structures of a non-rotating dipolar condensate. 

The density oscillations of the vortex state also af- 
fect its stability property. A condensate with sufficiently 
larg e and negative a is unstable and tends to collapse 
[i0|. The critical scattering length a cr as a function of 
the trap aspect ratio is shown in Fig. |2tb) for both the 
axial vortex state and the non- vortex state. Near a cr , 
due to the dipole-induced density oscillation, the vortex 
state has a peak density exceeds that of the non-vortex 
state. Since collapse starts locally at the high density 
region, the vortex state tends to be less stable, having a 
critical scattering length smaller in magnitude by about 
5~10% compared with the non- vortex state. This is in 
contrast with the non-dipolar condensate where the vor- 
tex states are shown to be more stable ^lj as they have 
lower peak densities compared to the non-vortex states. 
We remark that the repulsive dipolar interaction in this 
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FIG. 3: (Color online). Density profile of a vortex lattice state 
for A = 100, = 0Alu±, and a — — IOceb. The inset shows the 
position of the vortices which form a triangular lattice with 
hexagonal symmetry. 



case helps to stabilize a condensate with attractive scat- 
tering length: with all the other parameters being the 
same as in Fig. IHb), a non-dipolar condensate would 
have a critical scattering length on the order of — Chios . 

As we increase the rotation frequency, more vortices 
enter the condensate and form a regular lattice as a re- 
sult of the repulsive vortex-vortex interaction. For a 
non-dipolar condensate, it is well known that the vor- 
tex lattice takes a triangular shape with hexagonal sym- 
metry. Very recently, it is predicted by several authors 
that the vortex lattice of a dipolar condensate in the 
high- rota ting quantum Hall limit may take different ge- 
ometries ) 1 2L Il8j. However, we find from our numerical 
calculations that the triangular lattice with hexagonal 
symmetry always has the lowest energy for rotating fre- 
quencies up to Q — 0.99lu±. Although computing power 
and numerical accuracy currently prevent us from taking 
higher values of f2, this suggests that the quantum Hall 
regime of a fast rotating dipolar condensate may need a 
more careful analysis. An example of the vortex lattice is 
shown in Fig.[3]where the crater structure is still present. 
To ensure that such a structure is not an artefact due to 
the reduction of the effective spatial dimension, we also 
performed full 3D calculations and found that the results 
are in complete agreement with the 2D calculation. 

Transversely polarized dipoles — It is also of interest 
to study the vortex state with dipole moments polarized 
transversely. We assume that the dipoles are polarized by 
a transverse magnetic field co-rotating with the conden- 
sate about the z axis 0| . Without loss of generality, we 
assume that the dipoles are polarized along the x axis in 
the rotating frame. The effective quasi-2D dipolar inter- 
action potential for this case has a similar form as Eq. (0 
with f^(g) replaced by 
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FIG. 4: (Color online). Structure of a single vortex state 
with dipoles polarized along the x axis. The parameters are 
A = 100, a = 16as and Q. = 0.3uj±. (a) Density profile. The 
solid curves show the densities along the x- and j/-axis. (b) 
Density contour plot where the elliptical shaped vortex core 
can be clearly seen. Brighter color represents higher density, 
(c) Contour plot of the phase of the wave function. 



is anisotropic: it is repulsive along y axis and attractive 
along x axis. Such a system is unstable for scattering 
lengths below a threshold value, which is about 15as for 
the parameters used in our calculation. 

Figure0|shows the structure of a single vortex state for 
a = 16as. The whole atomic cloud in this case is elon- 
gated along the x axis with a two-fold symmetry. This 
is due to the magnetostriction induced by the dipolar in- 
teraction, an effect that has been recently observed in 
experiment ^3 . One can also notice from Fig. 0Jb) that 
the vortex core is also anisotropic: it has an elliptical 
shape with the major axis along x. This can be under- 
stood as follows: the vortex core size is determined by 
the healing length, and the attractive (repulsive) dipolar 
interaction along x (y) direction weakens (enhances) the 
contact interaction, resulting in a larger (smaller) effec- 
tive healing length, and hence core size, along x (y). 

Along any closed curves around the anisotropic vor- 
tex core, there is still a phase slip of 2tt. However, the 
phase of the wave function no longer coincides with the 
azimuthal angle ip [see Fig. EI C )]- Generally, a vortex 
state with a two-fold symmetry can be represented as 
if> ~ ^ n a n yo' 2 ™ +1 'e^ 2 " +1 ) v . To a good approximation 
the vortex core structure shown in Fig. 0] can be mod- 
elled by the three most dominant terms with n = 0, ±1, 
i.e., the wave function near the core takes the form 



It can be easily shown that, unlike in the previous 
case, now the effective 2D dipolar interaction potential 



4>{p) 



p(e llp + ae 



VP 



3 e 3iip 



4 




FIG. 5: (Color online). Density contour plots for A = 100, 
a = 18as and various rotation frequencies. The dots at the 
lower right corner shows the position of vortices. For SI = 0.5 
and 0.6ojx, the angles i) are 55° and 62°, respectively. 



with parameters a < and j3 > 0. 

At higher rotation frequencies, vortex lattice will form 
with a structure strongly depending on the scattering 
length and the rotating frequency. The richness of the 
vortex lattice structures are illustrated in Fig. [5] for 
a = 18ag at three different rotation frequencies. For 
SI = 0.35oj±, three vortices appear, all located along the 
x axis. Note that for the isotropic interaction as in the 
previous case or in a non-dipolar condensate, the three- 
vortex state would have the three vortices situated on the 
corners of an equilateral triangle. At higher rotation fre- 
quency, the vortices form lattices with two-fold symmetry 
instead of the usual hexagonal symmetry. For SI = 0.5 
and 0.6w_l, the angles ■& as defined in the middle plot of 
Fig. [5] are about 55° and 62° , respectively. Numerically, 
it becomes increasingly difficult to find the lowest energy 
lattice structure when we further increase SI. We plan to 
study this problem in greater detail in a future work. 

We remark that anisotropic vortices have been exten- 
sively studied in the context of high T c superconduc- 
tors, where the anisotropy can arise from the crystal 
anisotropy 16] or from high angular momentum pair- 
ing interaction (e.g., in d-wave superconductors) [l7| . In 
these anisotropic superconductors, the vortex lattice also 
deviates from the hexagonal geometry and can sometimes 
take a square shape jlq . 

In conclusion, we have investigated the properties of 
the vortex states of a quasi-2D rotating dipolar conden- 
sate. The vortex states possess many novel and unique 
features due to the presence of the dipolar interaction. 
We hope that our work will stimulate more experimental 
efforts on rotating dipolar quantum gases. 

This work is supported by Rice University and the Oak 
Ridge Associated Universities. 
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